Nonperturbative semiclassical stability of de Sitter spacetime for small metric 

deviations 
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We consider the linearized semiclassical Einstein equations for small deviations around de Sitter 
spacetime including the vacuum polarization effects of conformal fields. Employing the method of 
order reduction, we find the exact solutions for general metric perturbations (of scalar, vector and 
tensor type) . Our exact (nonperturbative) solutions show clearly that in this case de Sitter is stable 
with respect to small metric deviations and a late-time attractor. Furthermore, they also reveal a 
breakdown of perturbative solutions for a sufficiently long evolution inside the horizon. Our results 
are valid for any conformal theory, even self-interacting ones with arbitrarily strong coupling. 
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I. INTRODUCTION 

De Sitter space can be understood as an exponen- 
tially expanding cosmological spacetime entirely driven 
by a positive cosmological constant. It plays a central 
role in most models of cosmological inflation, where the 
potential-dominated energy density and pressure of the 
inflaton field act approximately as a cosmological con- 
stant and lead to a quasi-exponential accelerated expan- 
sion. The inflationary scenario provides a natural mecha- 
nism, through parametric amplification of quantum vac- 
uum fluctuations, for generating a nearly scale-invariant 
spectrum of adiabatic primordial inhomogeneities, which 
can successfully explain the observed CMB anisotropics 
and the large scale structure of the universe [l|, l2l • Fur- 
thermore, the physics of de Sitter could also be important 
for elucidating the final fate of the universe if its current 
accelerated expansion is entirely due to a small cosmolog- 
ical constant, a possibility compatible with observations 
so far. 

The exponential expansion quickly redshifts away any 
initial perturbations. This offers a simple means of estab- 
lishing natural initial conditions for subsequent evolution 
once such an accelerated expansion has already started 
over a region with a size larger than the Hubble radius 
and lasts for a sufficiently large number of e-foldings. Un- 
der these conditions the initial classical perturbations are 
effectively erased and the quantum state for modes with 
wavelengths much smaller than the Hubble radius is very 
close to the Bunch-Davies or Euclidean vacuum, which 
locally is essentially equivalent to the Minkowski vacuum 
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since at those length-scales the spacetime appears almost 
flat. Such a scenario and the late-time attractor charac- 
ter of local de Sitter spacetime, often referred to as the 
"no-hair" property of de Sitter, is supported by a num- 
ber of results and theorems in classical general relativity, 
both for linear perturbations [3, H| as well as for the full 
nonlinear case [5H8J. 

It is, therefore, of great interest to establish whether 
those classical no-hair results can be extended to the 
quantum mechanical case. Solid conclusions have re- 
cently been obtained within the framework of quantum 
field theory in curved spacetime 0, [1(3]. Specifically, 
given a fairly general class of massive interacting theo- 
ries with sufficiently weak coupling evolving on a fixed 
(nondynamical) dc Sitter background, it has been shown 
to all orders in perturbation theory [TJ [l2[ that quantum 
correlators within a spacetime region of bounded physical 
size become at sufficiently late times arbitrarily close to 
those of the Euclidean vacuum (the generalization of the 
de Sitter-invariant Bunch-Davies vacuum to interacting 
theories) . 

Considering test fields evolving on a fixed background, 
however, offers an incomplete answer: addressing the full 
dynamical problem requires taking into account the back- 
reaction of the quantum fields on the dynamics of the 
spacetime geometry. A number of studies have explored 
this question in the context of semiclassical gravity, where 
the metric is still treated classically, but its dynamics 
is governed by a generalization of the Einstein equation 
which includes the expectation value of the stress ten- 
sor operator of the quantum matter fields as a source 
[ToL [l3j . Focusing on the backreaction of conformal mat- 
ter fields, the dynamics of scalar-type metric perturba- 
tions around de Sitter has been analyzed in ref. flij . 
where the importance of considering also perturbations 
of the initial state of the matter fields has been empha- 
sized. (The linear stability of de Sitter including initial 
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classical stress tensor sources had earlier been studied in 
refs. [H|,[l6j].) In addition, the evolution of tensor metric 
perturbations for the same situation has been studied in 
ref. [13], where the semiclassical Einstein equation was 
solved perturbatively. (As we will show below, however, 
these perturbative solutions cease to be valid for a suffi- 
ciently long evolution inside the horizon.) The stability 
of tensor perturbations including the backreaction from 
conformal fields has also been considered [HI in investi- 
gations on the existence of unstable runaway solutions of 
the corresponding higher-order equations in the context 
of inflationary models driven by the trace anomaly [l9| . 
Nevertheless, those analysis neglected the contribution 
of nonlocal terms which play a key role in the existence 
of runaway solution for perturbations around flat space 

MM- 

Related studies have also been carried out for noncon- 
formal fields in de Sitter spacetime. The case of massless 
minimally coupled free scalar fields has been considered 
in ref. [21] and a vanishing correction to the classical 
modes was found when solving perturbatively the lin- 
earized semiclassical equation for tensor perturbations. 
Both massless and massive nonconformal scalar fields 
were studied in refs. [H, HH and the stability of de Sitter 
spacetime with respect to spatially isotropic perturba- 
tions was established. A fairly general class of Gaussian 
initial states was considered and de Sitter was found to be 
a late-time attractor in all cases. Moreover, the impor- 
tance of taking into account the contribution of nonlocal 
terms for light massive fields when analyzing the stability 
in the infrared regime was elucidated [23| . 

Here we consider the linearized semiclassical Einstein 
equation around a de Sitter background including the 
vacuum polarization effects of conformal matter fields, 
and solve it exactly for general metric perturbations (of 
scalar, vector and tensor type). In doing so, we make use 
of the method of order reduction [l3l . l24j , which elimi- 
nates the spurious solutions associated with higher-order 
derivatives while capturing the right dynamics in the in- 
frared regime. Moreover, the method generates a backre- 
action equation which is equivalent to the original semi- 
classical Einstein equation up to the same order in in- 
verse powers of the Planck mass at which the latter is 
valid within an effective field theory (EFT) approach to 
perturbative quantum gravity [25L [2q , but which can be 
significantly simpler to solve, a fact that we exploit in 
our calculation. Our exact solutions clearly show that 
de Sitter spacetime is also stable in this case and a late- 
time attractor as far as local geometrical properties are 
concerned. Furthermore, it reveals a breakdown of per- 
turbation theory when solving the semiclassical equation 
for a long time evolution inside the horizon. 

It should be stressed that our results are valid for any 
conformal field theory (CFT), even self-interacting ones 
with arbitrary strong coupling, as explained in sec. HXl In 
addition to metric perturbations, perturbed initial states 
of the matter fields have also been considered. 

Although semiclassical gravity does take into account 



the backreaction of the quantum matter fields on the dy- 
namics of the mean spacetime geometry, it does not pro- 
vide a complete analysis because it does not include the 
quantum mechanical effects of the metric itself. Indeed, 
one needs to quantize the metric perturbations in order 
to account for certain relevant phenomena: doing so is 
necessary, for instance, for a proper description of the 
generation of primordial cosmological perturbations, and 
it has even been suggested that radiative corrections in- 
volving higher-order graviton loops could lead to a secu- 
lar screening of the cosmological constant [13, [H| . How- 
ever, detailed calculations including graviton loops are 
technically complex and one is, in addition, confronted 
with the need to consider appropriate observables which 
are not only g aug e-invariant beyond linear order but also 
infrared safe [29143 lj . Because of such difficulties only 
partial progress has been made in this direction. It is, 
therefore, important to consider also somewhat less ambi- 
tious problems, but obtain solid results (and, if possible, 
exact) which can provide a robust foundation for further 
developments. One such example is the exact calculation 
of one-loop corrections from matter fields to the correla- 
tor of the Riemann tensor [32-34] for quantized metric 
perturbations around de Sitter. The results support the 
existence of quantum states for metric perturbations in- 
teracting with matter fields which exhibit (appropriately 
defined) de Sitter invariance, at least when graviton loops 
are neglected. In contrast, the results presented here 
only apply to the mean field geometry, but explore the 
effect of different (non-de Sitter-invariant) initial states 
on the dynamics and show not only that a self-consistent 
de Sitter-invariant solution exists, but also that it is a 
late-time attractor. Furthermore, the methods described 
below for obtaining nonperturbative solutions valid for 
long evolution times could prove helpful in order to ex- 
tend the calculation of the Riemann correlator so that 
it correctly captures the details of its behavior for large 
separations (both spatial and temporal), which seems to 
require a nonperturbative treatment. 

The rest of the paper is organized as follows. Semiclas- 
sical gravity and the method of order reduction are briefly 
reviewed in sec. [TTJ The semiclassical Einstein equation 
for metric perturbations around a spatially flat FLRW 
spacetime including the quantum back-reaction of con- 
formal fields is presented in sec. lIIII Given a cosmological 
constant and fields in the Bunch-Davies vacuum, there is 
a self-consistent semiclassical de Sitter background. In 
sec. IIVI the metric perturbations are decomposed into 
scalar, vector and tensor contributions. Next, we fix the 
gauge and write the (decoupled) semiclassical equations 
for the three types of perturbations, before and after em- 
ploying the order reduction method. The exact (non- 
perturbative) solutions are obtained in sec. [V] and their 
implications for the stability of de Sitter spacetime as 
well as the breakdown of the perturbative solutions are 
analyzed in sec. lVIl The effects of perturbing also the ini- 
tial state of the matter fields are studied in sec. IVIII and 
we show that all our main conclusions remain unchanged. 
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Finally, in sec. IIXI we summarize and discuss our results 
as well as explaining their applicability to any CFT for 
the matter fields. Several useful formulae concerning the 
perturbative expansion of curvature tensors, their con- 
formal transformations and some special functions are 
provided in the first three appendices. In addition, a 
method for generating a family of regular Gaussian ini- 
tial states is described in appendix [Dl and in appendix [E] 
we compare our linearized semiclassical Einstein equation 
for tensor perturbations (before order reduction) with the 
one previously obtained by Starobinsky [35j . 

Throughout the paper we use natural units with c = 
h = 1 and take k 2 = 167tGn- We employ the "+++" sign 
convention of ref. [3(| and Greek indices range over space 
and time, while Latin indices denote spatial components 
only. 



II. SEMICLASSICAL GRAVITY AND ORDER 
REDUCTION 

A. Semiclassical Einstein equations 

Semiclassical gravity can be regarded as a mean field 
approximation to a quantum theory of gravity where the 
mean gravitational field is treated classically and only 
matter is quantized. In contrast to quantum field theory 
in curved spacetime, it also includes the back-reaction 
of the matter fields on the mean geometry, given by 
the background metric g^ v . To achieve this, the stress 
tensor on the right-hand side of the Einstein equation 
is replaced by the expectation value of an appropriate 
quantum stress tensor operator. This expectation value 
needs to be renormalized, and counterterms local in the 
gravitational field have to be included in the bare grav- 
itational action. At all loop order for the matter fields 
(but no graviton loops) those counterterms are quadratic 
in the curvature and the renormalized semiclassical Ein- 
stein equation reads 

G^y + A(fi)g^ = ai(/i)A M „ + a 2 {n)B f _ ll , 

Z \ / ren 

where = (-g)~ 1/2 ( W") J C a ^ 5 C aMS ^d 4 x 
and = {-g)~ 1/2 (5/5g» v ) J R 2 y/=gd 4 x are obtained 
by functionally differentiating the finite parts of the grav- 
itational counterterms. The parameters ai(/x) and a2(p) 
together with A(/i) and k 2 (/j,) are in general renormal- 
ized parameters which have to be determined by exper- 
iment, and fi is the renormalization scale. Note, never- 
theless, that the backreaction equation ([1]) is renormal- 
ization group invariant and the dependence on \i of the 
different parameters appearing in the equation and the 
renormalized expectation value (T /1I/ (/i)) ren cancel out. 



The tensors A^ v and B^ are explicitly given by 
V - -4V( a V«C ai ^ - 2R^C a ^ 

— —4RanR" + ^RRfiv — ^g^R 2 + g^R" 13 R a p 
+ 8V [a V M]j R° - 2D g Rp, + ^V M V„i? + -g^D g R 
B^ v = -g^uR 2 — 2RR^ V + 2V M V„i? — 2g^D g R, 

(2) 

where the two forms of A^ v given above were obtained 
by using the identity (|A5I) and the definition of the Weyl 
tensor (|A4j) as well as the second Bianchi identity. 

B. Order reduction 

The semiclassical Einstein equation (JlJ contains terms 
with up to fourth-order derivatives of the metric, as seen 
from cq. ([2]). (The expectation value (T M „(/z)) ren also 
involves similar terms, as shown below.) Such kind of 
higher-order time derivatives are common in backreac- 
tion problems. A well known example is the Abraham- 
Lorentz-Dirac equation, which describes the effect of ra- 
diation reaction on the motion of a point-like charge in 
classical electrodynamics [33, HH (i.e. without consider- 
ing the internal structure of the particle nor a finite size 
for the charge density distribution). In fact, they are a 
generic feature of effective field theories (EFTs), where 
the effects of the UV sector on the dynamics of the low- 
energy degrees of freedom are encoded at the level of the 
action through an expansion of local terms with an in- 
creasing number of derivatives. The validity of the EFT 
expansion relies on the fact that for length-scales much 
larger than the inverse cut-off scale of the UV sector the 
higher-order terms in the expansion become increasingly 
smaller. In this regime their contribution amounts to a 
small correction to the equation of motion which results, 
when treated perturbatively, into locally small perturba- 
tions of the classical solutions. In contrast, solving the 
corresponding higher-order equations exactly gives rise to 
additional solutions exhibiting exponential instabilities 
with characteristic time-scales comparable to the inverse 
cutoff scale of the EFT (or sometimes fast oscillations 
with the same kind of characteristic timescale) , often re- 
ferred to as "runaway" solutions. These are spurious so- 
lutions which should not be taken seriously since they in- 
volve characteristic scales for which the EFT expansion 
breaks down and the contributions from the higher-order 
terms to the equation of motion no longer correspond to 
small corrections but to dominant terms. 

The simplest way of avoiding such spurious solutions 
is by solving the corrected equations of motion pertur- 
batively. However, perturbative solutions may not be 
valid for long times. This happens when quantities like 
the total time appear multiplying the perturbative pa- 
rameter so that the expansion contains so-called secular 
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terms which grow with time and lead to a breakdown for 
sufficiently long times of the truncated perturbative ex- 
pansion. Those limitations can be overcome with the or- 
der reduction method, which consists in taking the equa- 
tion of motion with corrections up to a finite order and 
writing an alternative equation which is equivalent up to 
that order but contains no higher derivative terms (this 
is achieved by taking successive derivatives of the original 
equation and substituting the higher-order derivatives in 
the correction terms to the appropriate order). The ex- 
act solutions of the equation obtained with this method 
agree locally with the perturbative solutions constructed 
around different times (each one with a finite domain 
of validity) and provides an interpolation between all of 
them valid for long times. This is particularly impor- 
tant when considering situations where the effects of the 
corrections are locally small, but can build up over long 
times and give rise to substantial accumulated effects. 
Two examples of such situations are an electric charge 
following a quasi-circular trajectory in a uniform mag- 
netic field and emitting electromagnetic radiation for a 
sufficiently long time so that the radius of its orbit de- 
creases, say, to half of its initial value due to radiation 
reaction, or an evaporating black hole emitting Hawking 
radiation for such a long time that its mass (or horizon 
size) decreases to a small fraction of its initial value. 

Related alternative methods which have been em- 
ployed in the literature for discarding the spurious solu- 
tions mentioned above involve finding the exact solutions 
of the original backreaction equation and then selecting 
the appropriate subset either by demanding analyticity of 
the solutions with respect to the perturbative parameter 
or checking explicitly which solutions exhibit unphysi- 
cal characteristic scales and disregarding them. How- 
ever, the latter method is less systematic and requires 
a case- by-case analysis, whereas the analyticity require- 
ment may be too restrictive in some cases (l3| . Further- 
more, the order reduction method leads to equations of 
motion which are equivalent up to the order under con- 
sideration, but are often easier to solve, as will be the 
case for the problem analyzed in the remaining sections. 

The order reduction method has been applied to elec- 
tromagnetic [HI and gravitational [39| radiation reaction 
problems as well as higher derivative gravity 14011 . It 
has also been employed in semiclassical gravity [l3l [24| 
and in this context it has been argued [4l| that trace- 
anomaly-driven inflationary models (with no cosmologi- 
cal constant and driven entirely by the vacuum polariza- 
tion of large number of matter fields [19]) correspond to 
spurious solutions which lie beyond the EFT's domain of 
applicability and are automatically discarded when using 
order reduction. 

It should be noted that order reduction cannot be al- 
ways applied in a straightforward way. It may be am- 
biguous in integro-differential equations, or may lead to 
covariance breaking if the time derivatives and spatial 
derivatives are not simultaneously reduced; see [l3j for a 



detailed discussion of these issues. 

The order reduction method can be illustrated in a 
nutshell with the following simple example of a first order 
differential equation in time for a function f(rj) with a 
perturbative correction of order k 2 . Given 



/' + bf = K 2 P(fJ'j" ,...), 



(3) 



where & is a constant and P is an arbitrary function, 
order reduction uses that /' = —bf + C(k 2 ), and by 
deriving one more time /" = —bf + C(k 2 ) = b 2 f + 
C(k 2 ). Substituting those two equations into the right 
hand side, we get 



f + bf = K 2 P(f, -bf, b 2 f, ...) + 0{K A ) 



(4) 



which is an equation of first order which is valid to the 
same order in k 2 as the original equation ([3]), but does 
not have unphysical solutions. Rather than considering a 
truncated perturbative expansion, this equation can now 
be solved exactly. It is clear how the method works for 
equations of more derivatives or partial differential equa- 
tions: one takes the lowest order equation and substitutes 
it in the higher order terms (in k 2 ), taking additional 
derivatives if necessary. 



III. CONFORMAL FIELDS IN A PERTURBED 
FLRW UNIVERSE 

A. The model 

Our model consists of N massless free scalar fields con- 
formally coupled to the spacetime curvature: 



S[<p,g] 



'd^tpj d v tpj + Z, cc R<p) 



gd d x, 
(5) 

where R is the Ricci scalar associated with the metric g^ 
and £ cc = (d— 2)/4(tZ — 1), which reduces to £ cc — 1/6 in 
four dimensions. We will use dimensional regularization, 
but after the renormalization procedure has been carried 
out we will take d = 4. Furthermore, we will specialize 
the physical metric g^ v to a slightly perturbed spatially 
flat FLRW spacetime, which is conformal to an almost 
flat metric g M „: 



(6) 



where r\ denotes the conformal time. 

The renormalized semiclassical Einstein equation ([1]) 
for this model with the fields in the conformal vac- 
uum state was derived by Campos and Verdaguer in 
refs. [HI, using the closed-time-path (CTP) effective 
action, and is given to linear order in the perturbation 
h^ v by 
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Gfj, v + Ag^ = --^(3k 2 B^ + ^an 2 - 2R al3 C fiai/0 



2—2 

-an a 



-4V Q V^ (C^up lna) + J H(x-y;fl) A^&^y 



(7) 



where a — iV/(28807r 2 ). Here and throughout the rest 
of the paper denotes the covariant derivative with re- 
spect to the metric g^ v and D g = V M V M , whereas objects 
with a tilde are evaluated with the conformally related 
metric <? M „. Furthermore, background quantities will be 

denoted by a superscript ^ such as gj^J = a 2 ^)^^, and 
quantities linearized in the perturbation by a superscript 
(1) as in g$ = a 2 (77)/^. 

For conformal fields the renormalized parameter ei2 (/i) 
does not depend on fi and we denote it by /?. Moreover we 
have chosen a renormalization scale \x such that a\ {fi) = 
0. Similarly, both A and k 2 are also independent of /1 
in this case. The tensors and -B M „ are defined in 
equation ^ and 



H, 



1 



1 



(8) 



Finally, the kernel H(x — y; fi) depending on the renor- 
malization scale [i is given in appendix [El but its ex- 
act form will not be needed in the bulk of the paper. 
Eq. ([71 coincides with those derived by alternative meth- 
ods SHEl. 



Hence, de Sitter spacetime, given here in spatially flat 
coordinates (the Poincare patch), is a self-consistent solu- 
tion of the semiclassical Friedmann equation (J9j) , with the 
effective cosmological constant (fTTj) having a small posi- 
tive shift of quantum origin. The existence of such self- 
consistent solutions follows straightforwardly from the 
fact that the renormalized expectation value of the stress 
tensor for the Bunch-Davics vacuum must be propor- 
tional to the metric, as implied by de Sitter invariance, 
and has been know for a long time [ill [4?]]. When A = 0, 
eq. © still admits a de Sitter solution with H 2 = 2/ (an 2 ) 
(closely connected to Starobinsky's original model of in- 
flation pjl, [4|| ) , but its characteristic scale lies beyond 
the domain of validity of semiclassical gravity when re- 
garded as part of an EFT approach to quantum gravity, 
as briefly discussed in sec. Ill Bl Such solutions are au- 
tomatically discarded by the method of order reduction 
241. 



IV. LINEAR PERTURBATIONS 



A. Gauge fixing 



B. The FLRW background 

The semiclassical generalization of the Friedmann 
equation can be obtained by setting the perturbation h^ v 
to zero in the 00 component of eq. ([7]), which gives 

6(a') 2 - 2Aa 4 = 3a K 2 a- 4 {a') 4 

+ 3f3 K 2 a~ 3 [2aa'a"' - a(a") 2 - 4(a') V' 1 



and using the order reduction method it becomes 



(a'f 



A 



1 + -an 2 A 
6 



0{n 4 ) 



Defining an effective cosmological constant A e ff as 

1 



Aeff = A 
Eq. (fTO)) has the solution 



1 + -an 2 A 
6 



(9) 



(10) 



(11) 



(12) 



where the Hubble parameter H is given by 3H 2 = A e jj, 
and — 00 < 1] < 0. This solution is unique up to a shift 
of the origin of conformal time, 77 — >■ 77 — 770, and its sign. 



When considering metric perturbations around a given 
background, there is a gauge freedom (corresponding to 
local diffcomorphisms) associated with the mapping be- 
tween the background and the perturbed geometry. In- 
finitessimal diffeomorphisms generated by an arbitrary 
vector field induce the following gauge transformation 
of the physical perturbation a 2 h^ v introduced in Eq. ([6]): 



i 2 h 



i 2 h 



(13) 



Rescaling the arbitrary vector field £ M = a 2 £ M and using 
Eq. (fT2"]) for the scale factor a(r\) 1 Eq. (flU]) becomes 

2 

h^u -> hfu, + <9^£„ + <9„£ M + ~ vfo ■ (14) 

Before proceeding any further, it is convenient to decom- 
pose the perturbation h^ v exploiting the fact that the 
spatial sections of the background metric are maximally 
symmetric spaces [49Tj5l1 |. The spatial part hij decom- 
poses as 

= hj? + 2<9 (i wJ ) + didja + tS^ , (15) 

where 6' •<),!, ), v = = fthj? and 6"0,ir [ , = 0. The 
temporal components can be similarly decomposed as 

hoi = vj + dit/j , h 00 = (f> , (16) 
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where S^divJ — 0. In total, we have four scalars <fi, 
ip, a, t, two transverse vectors wj and vj (with two 
independent components each) and a transverse traceless 
tensor hj^ (with two independent components as well). 
Decomposing also the spatial part of the vector field £ M 
as 



6 



(17) 



where S^di^J 



0, wc can sec the behavior of the various 



components under a gauge transformation: 



, tt 

Hj 

a 

„.T 



'ho 

T 

Vi - 



til 



r -> r + -Co 
»7 

^ -> *A + £' + Co 



(18) 



2 ^o - -Co 



where primes denotes derivatives with respect to the con- 
formal time r\. Choosing £f, £ and £ appropriately, we 
can set wj , a and r to zero. We will work in this gauge, 
where the perturbation of the spatial metric is entirely 
given by the tensorial component: 



(19) 



This is the transverse traceless gauge, also known as spa- 
tially flat gauge when focusing on the scalar perturba- 
tions. This fixes completely the gauge if we restrict our- 
selves to metric perturbations that fall off at spatial in- 
finity. 



If we consider perturbations which do not necessarily 
fall off, there is still some residual gauge freedom which 
is not fixed by condition (|19[) . On one hand, there are 
transformations which leave hij invariant. One possibil- 
ity are those generated by £j which are functions only of 
the conformal time; this corresponds to translations on 
the spatial sections which can change for each surface of 
the foliation and lead to changes of vj . A second possibil- 
ity involves transformations generated by £j = b{rf)SijX J , 
£o = —b(rj)r)', this corresponds to dilations on the spatial 
sections while changing at the same time the surface of 
the foliation so that the expansion of the FLRW back- 
ground compensates for that, and leads to changes of <f) 
and ip while leaving hij invariant. On the other hand, 
the transformations generated by £j — EijX^ with in- 
constant and traceless induce changes of hij but leave it 
transverse and traceless. These residual gauge transfor- 
mations will play a role in sec. [V] to show that certain 
solutions are pure gauge. 



B. Semiclassical equations of motion 

Using the metric decomposition and the gauge fix- 
ing introduced in the previous subsection, one can ob- 
tain from the semiclassical equation Q the dynamical 
equations for the scalar, vector and tensor perturbations, 
which decouple form each other. From the 00 component 
of eq. (0 and the scalar part of its Oi component, one gets 
the following two equations for the scalar perturbations 
i\) and 4>: 



J 



0= -(3-5an 2 H z )H 2 a 



2„4, 



2 [1 - (3a - P) n 2 H 2 } Ha 3 A^j + i/3 K 2 



7H 2 a 2 A6 + 3H 2 a 



2„2j« 



6ii-V 



3„3 il 



-A 



/3k 2 



- HaAip'' 



HaA 2 ^ + -AV 



2A 2 V>V) - A 2 (j){x') 
= 2Ha 3 Acj) + 2an 2 (HaA 2 cj) - 3H 3 a 3 Acj) ~ 2HaA 2 ^') + (3k 2 



(H(x - x'; ft) + S 4 (x - x') In a) dV, 



HaUAcf) - 2H 2 a 2 A(t)' - -A 2 (j}' 

3 



/3k 2 ( -4H 2 a 2 A 2 tp + t^A V 



2an x 



2AVV) " AV (x') 



(H{x - x'; ft) + <5 4 (x - x') In a) dV, 



(20) 

where □ = rf" u d^d v and A — S^didj. Similarly, from the transverse part of the Oi component one gets the equation 
for the vector perturbation vj: 







[1 - („ + 2ffl - 3«^«A„f + W / CA« 7 (/,) (H(x - « fl ) + - z> .) dV. 



(21) 



Finally, the equation for the tensor perturbations is obtained from the transverse and traceless part of the ij compo- 
nents: 



= -2 [l - (a + 2/3) k 2 H 2 } Ha 3 ti^ T + [l - (a + 2/3) k 2 H 2 ] a 2 DhJ^ + 3an 2 H 2 a 2 [2h'lJ T + OhJ^] 
- 6an 2 HaOti^ T + 3an 2 [ (DDhJ^(x')) (H(x - x'\ ft) + 5 A (x - x') In a) dV. 



(22) 



Employing order reduction as explained in sec. IIIB1 
the equations can be rewritten in the much simpler form 
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/4? T -Z(l-v) h£ T (1 - 2.) = 



(23a) 



(23b) 



A</> = e>( K 4 ) 



(23c) 
(23d) 

where we have introduced the following parameter, which 
controls the expansion in powers of k 2 : 

v = 3an 2 H 2 < 1 . (24) 

It is worth emphasizing that those equations are inde- 
pendent of the arbitrary parameter j3 and the renormal- 
ization scale p, of the semiclassical theory. They involve 
only the semiclassical parameter a, which depends on the 
matter field content. 



Nonlocal terms 



so that g^A^ = 0. 

Expressing now the Weyl tensor in terms of the Rie- 
mann tensor and its contractions according to its defining 
equation (|A4I) , we can write 



F = 2 



Jtt — —H 



gd x + / Siyj—gd x. 



(29) 

where £4 is the integrand of the four-dimensional Euler 
invariant and is given by 



£4 — RWP° R — kR^R + R 2 



(30) 



The generalized Gaufi-Bonnet theorem establishes that 
the integral of £4 is a topological invariant, namely 327r 2 
times the Euler characteristic. From eq. ([29)) and the fact 
that the variational derivative of a topological invariant 
vanishes, we see that can be expressed entirely in 
terms of the Ricci tensor, the Ricci scalar and covariant 
derivatives acting on them. 

We are now ready to use order reduction. Taking into 
account that 



F,uv 



g^R + Ag^ = 0(k 2 ) 
R-AA = 0(k 2 ) 



(31) 



As we have calculated explicitly, when we use order re- 
dution the nonlocal terms do not contribute to the semi- 
classical equations of motion (|23l) for the perturbation 
h^w We now want to show that one can see this in 
general, without choosing a gauge or expanding the semi- 
classical equations explicitly in terms of the perturbation 
hfj, u . Basically this amounts to showing that A^, given 
by eq. is of order k 2 when using order reduction. In 
order to do so, it is convenient to consider its definition 
as a functional derivative of the integral of the square of 
the Weyl tensor 



1 SF 
~^g Sgv 



(25) 



with 



F 



C a ^ s C a ^ s ^d 4 x = / C a ^ 5 C, 



^a^sy —gd x , 

. ! - ( : ; 

where the last equality follows from the conformal invari- 
ance of the Weyl tensor with one raised index. This also 
implies that 

a 2 6F 



1 SF 
^9 Sg^ 



a z A„ 



(27) 



and we can equivalently show that A M „ is of order n 2 . 
Another consequence of the invariance of F under con- 
formal transformations of g^ v is the vanishing trace of 
A ni , '. 



Q _SF _ SF d]T _ 2 
Sa Sg^ v da 



gA^gT , (28) 



and substituting i? M „ = Ag^ into the expression for A^ v , 
the result can only be proportional to g^ v up to order 
k 2 . However, since g^A^ — 0, we conclude that A^ 
is of order k 2 when order reduction is employed. One 
can alternatively check this fact by substituting eq. (|3"Tj) 
into the explicit expression for A^ v in terms of the Ricci 
tensor in eq. ©• 

In conclusion, we see that we only have to consider the 
local terms in eq. ([7]) when using order reduction. 



NONPERTURBATIVE SOLUTIONS 



A. Scalar and vector perturbations 



The solutions of eq. (|23bl) for the components of the 
vector perturbation vj are arbitrary functions of time, 
which can be eliminated by a gauge transformation. In- 
deed, by using the residual gauge freedom described at 
the end of sec. IIV Al and choosing £j as an appropriate 
function of time only, we can set vj = 0. 

For (j>, the solution of eq. (|23c[) is also an arbitrary 
function of time. The solution for ip is then given by 



V> = /(*) 



1 

4// 



\+-v)^{t)v 2 



(32) 



Since ip enters into the perturbation h^ u only through 
a spatial derivative according to eq. (|16p . the arbitrary 
function f(t) does not change the perturbation and 
we can set it to zero. If we want to start with bounded 
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initial perturbations, we must exclude solutions which 
are unbounded and have to take = 0. 

On the other hand, if we had not excluded such un- 
bounded solutions, we would first have to choose £o ap- 
propriately to make <f> vanish, and would need to take 
a similar unbounded function £ = — 1/ (2rj)^QY 2 so that 
the combination didja + rSij which enters into the de- 
composition of the perturbation (|15l) still vanishes. The 
solution for ip is then an arbitrary function of time which 
we can set to zero as above. 

Thus, we see that when order reduction is employed, 
both vector and scalar parts after solving the constraints 
are pure gauge and can be eliminated by a residual 
gauge transformation of the kind mentioned at the end 
of sec. ITVAl 

Note that Anderson et al. [14J , who investigated scalar 
perturbations without using order reduction, also con- 
cluded that those perturbations (which they refer to as 
perturbations of the first kind) have to vanish because 
the corresponding solutions that they found lie outside 
the range of validity of the semiclassical theory. The 
gauge-invariant variables that they introduce, however, 
take a simple form in a gauge rather different from the 
one we employ, and so their intermediate expressions are 
not directly comparable to ours. 



B. Tensor perturbations 

To find the solutions for the tensor perturbations, we 
take the Fourier transform with respect to the spatial 
coordinates, 



[ 4(p).9^p) eipx 

s=± J 



d 3 p 

(2tt) 3 



(33) 



where e^-(p) are a pair of transverse and traceless ten- 
sors corresponding to two different polarizations. Equa- 
tion ()23a|) then becomes 



n 



)9± 



(l-2v)p 2 g ± =0( K 4 



(34) 



Setting w 2 = (1 — 2i>)p 2 , s — —uj-q and g± = S2 _I/ /±(s), 
this reduces to a Bessel equation for f±, whose general 
solution is 

9± = {-unf-" [Cfj f _„(-w»j) + CfY^i-ur,) 

'(35) 

where Cf and Cf are integration constants. 

For the particular case p = 0, which corresponds to 
no spatial dependence in position space, eq. (|23[) can be 
solved directly and the general solution is given by 



Dij(-ri) 6 - JV + E, 



(36) 



where Dij and Eij are traceless tensors with respect to 
the induced background metric of the spatially flat sec- 
tions, and independent of the spatial coordinates and the 



conformal time0 The first term on the right-hand side of 
eq. (|36|) corresponds to a Bianchi I anisotropic deforma- 
tion of de Sitter, whereas the second one is pure gauge 
and can be eliminated by the residual gauge transforma- 
tion generated by the transverse vector = JSya^' /2. 



VI. DE SITTER STABILITY AND SECULAR 
TERMS 

A. Stability with respect to linear perturbations 

In this section we analyze the stability of the semiclas- 
sical de Sitter geometry with respect to small metric per- 
turbations including the back-reaction due to quantum 
vacuum effects from conformal fields. We do so by fo- 
cusing on the evolution of the Riemann tensor associated 
with the linearly perturbed metric g^ v , which has a num- 
ber of appealing properties. First of all, the linear per- 
turbation around de Sitter of the Riemann tensor R a P ^$ 
with appropriately raised indices is a gauge- invariant ob- 
ject. This follows from the fact that for the unperturbed 
background it can be written as R^ al3 -ys = ^H 2 8^Sgj, 
whose Lie derivative with respect to an arbitrary vector 
field vanishes. Furthermore, with this index structure 
the components remain unchanged when rescaling by the 
same constant the basis vectors of the tangent space at a 
given point. This implies that the components coincide 
with those in the physical basis of orthonormal vectors 
{a~ 1 do, a~ 1 di} of the background metric. Finally, the 
Riemann tensor provides a suitable characterization of 
the local geometry, in terms of which the stability and 
attractor nature of semiclassical de Sitter spacetime can 
be naturally formulated, as further discussed at the end 
of this subsection. 

In terms of the metric perturbations the linearized Rie- 
mann tensor is given by 



= 2H 2 8f^ s] h m + 2ff 2 7 7 V [ V 1 ^9 [7 V 
+ 2H 2 r,8^r,^ (2d (M K )0 - h'^) . 

(37) 

Using the gauge transformation (|14p. one can explicitly 
check that it is indeed gauge invariant at linear order. 
Fourier transforming with respect to the spatial coordi- 
nates and specializing to tensor perturbations, we get 

■ q [a0] q [a/3] \ ipx d3 P 



R (1)af) lS = 2H 2 



(2tt) 3 



(38) 



1 These solutions can also be obtained by taking the limit u — > 
of eq. (135 I t after rewriting C\ = Ci/cj 3-2 " so that a finite 
non-vanishing limit is obtained for the solutions associated with 
the two integration constants. In addition, the different ways of 
taking the limit p — > of (p) and the possibility of considering 
arbitrary linear combinations gives rise to the general traceless 
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where 



S°Jok = -(e ± ) j k ri9'± 



(39) 



stki=P i P k(e ± )iv 2 g± + sis°Jc 



Hence, we can see that all the Riemann components can 
be written in terms of g± and g'± . Since everything that 
will be said is entirely equivalent for both polarizations, 
for ease of notation we will omit in the remainder of this 
section the subindices ± labeling the two transverse po- 
larizations associated with each momentum p. 

Let us consider first the evolution of the Riemann ten- 
sor for modes well outside the horizon, i.e. with \ujr]\ <C 1. 
In this case one needs to evaluate the Bessel functions in 
eq. ([35]) using eqs. (|C3j) and (|C5|) . which leads to 



+ G{urj) = const. + 0{ujti) , 

+ 0(ojii) = const. + Q{urj) . 

(40) 



-HI-"'- 

4r(i-„,2,-» 



Substituting into eq. (|39p we see that the components of 
the Riemann perturbation in a physical basis decay like 
l/a = —Hr\ or higher order at late times, i.e. in the limit 
7/ 0. 

On the other hand, for modes inside the horizon, with 
\(jJT]\ ^> 1, one can use eqs. (|C4[) and (|C6I) to see that g 
and g' are of the form 



7T 



(41) 



.g'-^-Hl'-qi + oiiH] 



times an oscillatory factor corresponding to a linear com- 
bination of sin(w?y) and cos(w?7). Thus, from eq. (I39[) it 
follows that inside the horizon the components of the 
Riemann perturbation oscillate with an amplitude that 
decays like l/a 1- ". 

Putting these results together we can conclude that 
de Sitter spacetime remains stable with respect to small 
metric perturbations of the semiclassical mean geome- 
try when the quantum back-reaction of conformal mat- 
ter fields is included. This is guaranteed by the fact 
that for any Fourier mode with comoving momentum 
p the perturbation of the Riemann tensor decays like 
1 /a u (times an oscillatory factor) when the correspond- 
ing physical wavelength 27ra/|p| is smaller than the de 
Sitter radius 1/H and like l/a when it is larger, together 
with the regularity of the perturbation around horizon 
crossing (when the wavelength is comparable to 1/H). 
This extends the conclusions of the no-hair theorem for 
de Sitter spacetime, which is not only stable with re- 
spect to small metric perturbations but also a late-time 
attractor in classical general relativity, to the case where 



radiative corrections from loops of conformal fields are 
considered. In fact, the main effect of the radiative cor- 
rections compared to the classical case for pure gravity, 
which corresponds to taking v — in our results, is sim- 
ply to alter slightly the exponent of the power-law decay 
for modes inside the horizon. 

Our result can be used to illustrate in a simple way 
the fact that the stability and the character of late-time 
attractor of de Sitter spacetime applies to sufficiently lo- 
calized observables characterizing the geometry within 
a region of fixed physical size (as opposed to comov- 
ing). The tensor perturbation /i^ T (?7,x) and the ampli- 
tude g(j],p) associated with a given momentum and po- 
larization are gauge-invariant objects with well-defined 
geometrical meaning. However, as their characteristic 
physical wavelength gets exponentially redshifted, at late 
times one would need to measure them over regions with 
a physical size that becomes arbitrarily large. Instead, 
the deviations of the geometric properties within a re- 
gion of fixed physical size compared to those of de Sitter 
decay exponentially with the (proper) cosmological time. 
These features are adequately captured by the behavior 
of the Riemann tensor, which provides a good character- 
ization of the local geometry. 



Perturbative vs. nonperturbative solutions 



In this subsection we will compare the exact nonper- 
turbative solutions of the linearized semiclassical equa- 
tion (|3"4"]l . given by eq. (|3"5"|) . to those that result from 
solving the equation perturbatively in v. We will see 
that the perturbative solutions cease to be valid when 
the modes evolve inside the horizon for a sufficiently long 
time. 

Here we concentrate on initial perturbations corre- 
sponding to Bunch-Davies positive frequency modes, but 
our conclusions can be easily generalized to arbitrary ini- 
tial conditions. This choice amounts to setting C2 = iCi, 
so that the linear combination within the square bracket 
on the right-hand side of eq. (|35p becomes a Hankel func- 
tion of the first kind, with a purely positive frequency os- 
cillatory behavior at early times. Indeed, using eqs. (IC4[) 
and (|C6|) one can see that in this case the behavior for 
\uj7]\ ^ 1 of the non-perturbative solution is given by 



(/ = r lA /l e - i t( 2 - 1 ') (-^'-"e-^ii + ofiH] . 

(42) 

On the other hand, if we first expand the solution ([33)1 
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in powers of v employing eqs. (|C8I) - (|C9I) . we get 



9 P = 



-Ci j|i (1 + i|p|r?) e-'IPl" - v (2 + pV) e^l" 
-\v [sin (|p|7 7 ) — |p|77Cos(|p|ry)] [in (p V) — i?r] 
- v (1 - i|p|?7) [Ein (-2i|p|?7) + 7 + In 2] e i|p| " 

-o(, 2 ). 

(43) 



At late times (for <C 1) the exact solution is well 
approximated by the perturbative solution (|43[) : one re- 
covers the asymptotic behavior in eq. f)40[) but with the 
constant term and the coefficients of the higher-order 
ones given by their expansion in powers of v truncated 
at linear order. That is, however, not the case at early 
times, with |w?7| 3> 1. This can be seen by considering 
the early-time limit of the perturbative solution (|43l) : 



9p 



\P\Vs 



1-y +ln(-|p|?/)-i|p|r? 



(44) 



which coincides with the result that one obtains by ex- 
panding in powers of v the asymptotic expression (|42[) 
for the exact solution. It is clear that the perturbative 
solution deviates significantly from the exact one when 
HpIM ^ 1- The implications can be more easily un- 
derstood if we normalize the modes so that they have 
a fixed amplitude at the initial time 770 independently 
of the particular value of 770 , which can be implemented 
by dividing the mode by its value at tjq. Proceeding in 
this way with the exact solution, one finds for example 
that the amplitude of a mode which was initially well 
within the horizon has decreased by a factor l/(— a;ryo) 1_1/ 
by the time of horizon crossing. In contrast, repeating 
the procedure with the perturbative solution and treat- 
ing v perturbatively when normalizing by the value at 

I 



7/0, one obtains an amplitude at horizon crossing of or- 
der (— l/|p|?7o) (l + i/ln(— |p|?7o) — HpI^o)- For modes 
which have spent a long time inside the horizon, so that 
Hp 11% I <^ 1j this amplitude at horizon crossing can be 
significantly larger. 

The reason for the potentially large deviation of the 
perturbative solution with respect to the exact one is a 
breakdown of perturbation theory: the actual condition 
for the validity of the perturbative solution is ^|p|?7o <C 1, 
which can be violated even for v <C 1 when consid- 
ering I p> ) J 770 1 large enough. The term proportional to 
if|p|?7 in eq. (|44l) . responsible for the main deviations, 
is a secular term arising from the truncated perturba- 
tive expansion in powers of v of the oscillatory factor 
e -iwr) j n ^ (j42|)_ Such a breakdown of perturbation 
theory for long times is very common when determin- 
ing the evolution of a system by solving perturbatively 
the corresponding dynamical equations. This can be 
illustrated with the simple example of a harmonic os- 
cillator with frequency + SQ. If one solves pertur- 
batively in SQ the corresponding equation of motion, 
x+(£l+6n) 2 x — 0, to first order one finds a solution of the 
form x(t) « Asm(flt + <p)[l + 6 fit], where the correction 
grows with time. In this case, however, the exact solu- 
tion is obviously known: x(t) — Asm [(SI + 6SY)t + <p\. It 
is, therefore, clear that the exact solution is qualitatively 
very similar to the unperturbed one but with a slightly 
corrected frequency. The growing terms in the perturba- 
tive solution, which are commonly know as secular terms, 
reflect the fact that the perturbative solution is not valid 
for arbitrarily long times but restricted instead to times 
such that Sflt <C 1. This is a rather simple example, but 
the situation is completely analogous for our semiclassi- 
cal solution. 

We close this subsection by comparing the perturba- 
tive solution obtained by expanding the exact solution, 
which has been discussed above, with the one obtained 
by solving perturbatively eq. Q34[) as done in ref. [17]. 
Substituting the classical solution for g into the terms 
proportional to v in eq. (|34[) . treating them as a source, 
and using the retarded propagator for the unperturbed 
equation (with v = 0), one gets the following result for 
the perturbative correction, already obtained in ref. fij\ 
(see their eq. (48)): 



cx 2i ^ 



[(l + P 2 77t/) sin [|p|(77 - 77')] - |p|(r? - 77') cos [\p\(rj - 7/)]] 



a-ilPlV 



■drj' . 



(45) 



The integral can be done exactly with the help of ap- 
pendix [Cj For large negative values of 770 one gets 



-i\p\ Vo + In (-|P|%) + |c- 2i IPl^ + 0(1) 



(46) 

where we have set 77 = for the upper limit, whose ef- 



fect is unimportant for a late-time solution. (Note that 
the logarithmic term and the phase were not included in 
ref. [17j . where only the dominant term was retained.) 
If we had normalized the classical solution at the initial 
time, which essentially amounts to dividing the unper- 
turbed solution by 1/tjq, the perturbative correction in 
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eq. (l46l) would agree with the perturbative solution (|44"]) 
when normalized at the initial time 770 as discussed below 
that equation. There is actually a slight discrepancy: one 
does not get the term e~ 2l l p l' ,0 /2 appearing in eq. ([4T>|). 
This is due to a different choice of initial conditions. 
Whereas our exact solution (|42j) results from choosing 
C2 = iCi exactly (at all orders in perturbation theory), 
eqs. (l4"51) - (|4l)f correspond to making the same choice for 
the zeroth-order solution but requiring the first-order cor- 
rection and its derivative to vanish at the initial time. By 
imposing the same conditions when determining our ex- 
act and perturbative solutions, one obtains a new pair of 
constants C\ and C2 which differ at order v and lead to 
a perturbative solution in full agreement with eq. (|46|) . 

It should be pointed out that the perturbative correc- 
tion can be large only when the corresponding mode was 
trans-Planckian at the initial time, i.e. its physical wave- 
length was much smaller than the Planck length (other- 
wise one has — f|p|?7o *C 1 and the perturbative correc- 
tion is always small), as recognized in ref. [l7|- At such 
scales semiclassical gravity is not guaranteed to provide 
an accurate description. Furthermore, one would need a 
rather small amplitude of the initial perturbations so that 
nonlinear gravitational effects do not become important: 
otherwise for such short wavelengths the effective stress 
tensor quadratic in the metric perturbations I53I ] can 
generate a strong back-reaction on the background ex- 
pansion, and even dominate over the cosmological con- 
stant. In any case, even if one carries out a linearized 

I 



where the the right-hand side corresponds to the pertur- 
bation of the expectation value of the stress tensor eval- 
uated on the background metric due to the perturbation 
of the initial state. Here it has been assumed that such a 
term is of the same order as the remaining terms, which 
are linear in the metric perturbations. Hence, no further 
terms should be considered since those corresponding to 
the perturbation of the initial state in the stress ten- 
sor expectation value evaluated on the perturbed metric 
would be of higher order. 

Note that eq. (|4"T1) and the procedure employed below 
can also be used to consider the effect of any other ad- 
ditional stress tensor sources (even classical ones) which 
can be treated perturbatively and regarded of the same 
order as the terms linear in the metric perturbations. 

For conformal fields in a FLRW background eq. (|47|) for 
the linear metric perturbations reduces again to eq. ([7]) 
plus the source term involving 8T^ W . Furthermore, one 
can introduce a decomposition of ST^ analogous to that 
for the metric perturbations in eqs. (j 15[) - f| 16[) . which 
is complemented now by the conservation requirement 
for STfw with respect to the background metric, i.e. 



analysis without much concern for these issues, as done 
in ref. [Tjj, the corrections to the classical solution will 
always be small as shown by our exact solutions and dis- 
cussed above. 



VII. PERTURBED INITIAL STATE 

So far we have considered the dynamics of metric per- 
turbations, but for a fixed initial state of the matter field, 
namely the Bunch-Davies vacuum. The state-dependent 
expectation value of the stress tensor is affected by the 
metric perturbations, but this is due to the their effect on 
the evolution of the scalar field operator in the Heisen- 
berg picture or, alternatively, on the evolution of the 
state in the Schrodinger picture. However, if one wants 
to allow changes in the initial state, eq. ([7]) needs to be 
generalized. In fact, as discussed in appendix [Dl the 
Bunch-Davies vacuum of the matter fields is no longer a 
Hadamard state (free of excitations at arbitrarily short 
wavelengths) for nonvanishing metric perturbations at 
the initial time: Therefore, the initial state needs to be 
modified so that it is a Hadamard state in that case and, 
in particular, the renormalized expectation value of the 
stress tensor is finite at the initial time. 

Starting with eq. (|TJ) and considering not only linear 
perturbations around the de Sitter metric but also small 
perturbations of the initial state, one obtains 



(47) 

I 

(ojy/Mjy^ = q w ith f'Vj, being the covariant deriva- 
tive associated with the background metric gjS) = a 2 ^^. 
More specifically, we can write 

ST 0l = 6T£ + d iX 

ST iS = 5T^ T + 2d (l W^ + Udj - I a) p + , 

(48) 

and take as independent quantities the transverse and 
traceless tensor ST? T , the transverse vector <5T(J and the 
two scalar functions <5Too and 5T. The temporal compo- 
nent of the conservation equation determines the scalar 
X and, similarly, the scalar part of the spatial projection 
of the conservation equation, which can be written as 
the gradient of a function, determines the scalar p. On 
the other hand, the transverse part of this spatial pro- 
jection determines the transverse vector Wj . Applying 
this decomposition to the generalization of eq. ([7]) in- 
cluding the source ST^, one obtains eqs. plus 
the corresponding sources. Finally, using order reduction 
eqs. (f23|) are generalized to 
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^ TT h'^ T - (1 - 2v) AhJ? = K *ST§* + O(^) , (49a) 

Ay? = -K 2 5TT i + 0(K 4 ) , (49b) 



A0 = -n 2 ( <5T 00 - ~?7<^ + ) + G(« 4 ) , (49c) 



A^+^ (1 + ^)0= ^« 2 <5T 00 + O(« 4 ) , (4!M) 



which correspond, respectively, to the transverse and 
traceless part of the ij component of the order-reduced 
equations, the transverse part of the Oi component, the 
spatial divergence of the Oi component, and the 00 com- 
ponent. Moreover, the conservation equation has been 
used to express the spatial divergence of ST oi in terms of 
STqo and ST on the right-hand side of eq. (|49cl) . 

A. Tensor perturbations 

Let us start with eq. (|49al) for the tensor perturbations. 
In addition to employing eq. (|33p for hj? , one can use 
the analogous expression 

«^T(ry,x) =J2j e-,(p)J^p)e ipx ^ (50) 

for the stress tensor perturbation and obtains the follow- 
ing equation for each one of the two polarizations: 

g'l - - (1 - u) g' ± + (1 - 2v) p 2 g ± = k 2 J ± . (51) 

This equation is identical to eq. (|34p but with an inhomo- 
geneous source term. The general solution can be written 
as g± = <j>k + gf , a sum of a homogeneous solution con- 
taining the information on the initial conditions at the 
time ?7o and a particular solution of the inhomogeneous 
equation with vanishing initial conditions, 

9Hv)= I G iet (r,, V ')J±(v')dv', (52) 
J no 

where G re t is the retarded propagator associated with 
eq. (|51[) . Given two independent homogeneous solutions, 
Mi and U2, the retarded propagator for such a linear 
second-order differential equation can be expressed as 

G Ist (r],r]) = 0(V-V), (53) 

where W{f}') — u[(r]')u2(r)') — u' 2 { , q')ui('q l ) is the Wron- 
skian for this pair of solutions, which is nonzero for in- 
dependent solutions. Therefore, the inhomogeneous so- 
lution can be written as 

9?(V) = Cf{v)uM - Cf{r,)u 2 {rf) (54) 



with 

c - M 'Cim JMW - (55) 

Provided that C^(rj) are regular and have a finite limit 
when rj — ¥ 0, as we will discuss next, the late-time behav- 
ior of the contribution form the inhomogeneous solution 
g^ is the same as for the homogeneous solution, already 
analyzed in sec. IVI Al We can consider the same pair of 
homogenous solutions as in that section and choose 

u 1 (i 1 ) = {-u:ri)i- v J?,_ v (-ujri), 

3_ (56) 
02(57) = (-W77) 2 "Y3_ u (-ur]) . 

Taking into account that for this choice W{rf) — 
(2lo/it)(— u-q') 2 ~ 2 ' J , we can immediately see that C^{rj) 
will be regular and have a finite limit when 77 — > pro- 
vided that J±{rj r ) is also regular and decays faster than 
(— r\') 1 ~ 2v as that limit is approached. This means that 
the components of the stress tensor ST^ T should decay 
slightly faster than 1/a in conformal coordinates, or 1/a 3 
in a physical basis. Such kind of behavior is fulfilled by 
classical radiation, which decays like 1/a 4 , and it seems 
plausible that small excitations of the vacuum state for 
the conformal fields considered here also exhibit the same 
decay; it is indeed the case for the class of regular states 
considered in appendix Q21 as shown there. We stress 
again that 5T^ V has been assumed to be small enough so 
that it is at most of the same order as the terms linear in 
the metric perturbation in eq. (|47l) and nonlinear grav- 
itational effects arc not important (otherwise one could 
consider for instance an excitation with sufficiently large 
energy overdensity to form a black hole by gravitational 
collapse, contradicting the conclusions above). 

As seen above, if the quantities in eq. ([5"5j) have a reg- 
ular limit as r\ — > 0, the inhomogeneous solutions gf^ 
behave in the same way as the homogeneous ones and 
the same conclusions drawn in sec. I VII apply to the inho- 
mogeneous solution as well. 

B. Scalar and vector perturbations 

Unlike tensor perturbations, vector and scalar ones 
do not give any nontrivial contribution in the classical 
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case, but they get a nonvanishing semiclassical correc- 
tion. Nevertheless, since they are governed by dynami- 
cal constraints which become simple algebraic equations 
when working in Fourier space for the spatial coordinates, 
their long time behavior can be directly determined from 
the fall-off properties of the stress tensor perturbation 

Let us consider the vector perturbations first. From 
eq. (|49bp it follows that in spatial Fourier space vj = 
(« 2 /p 2 ) <5Tq^ and the metric perturbations decay in the 
same way as the stress tensor. From eq. (|3T[) it follows 
that their contribution to the linear perturbation of the 
Riemann tensor is given by 



(57) 



and we can conclude that at late-times the curvature per- 
turbations decay like 1/aST^ and, hence, for any decay- 
ing (or even asymptotically constant) stress tensor per- 
turbation the stability and attractor character of de Sit- 
ter spacetime are not altered. 

The situation is similar for scalar perturbations but 
with some slight differences. From eqs. (|49c[) — ()49d|) one 
can see that </> decays at least like ST o or rj^ST^, 
whereas ip is only guaranteed to behave at late times like 
a5Too or aif^ST^. Given the contributions of scalar 
perturbations to the linearized Riemann tensor, 



H 2 r)p>pi (ip -r)ip') , 



R^ 0j k l=iH 2 ri5i k p l] <f > ., 



(58) 



one can see that they basically fall off like ST^v at late 
times. Therefore, if ST^ decays at least like l/o the 
conclusions about the stability of de Sitter space remain 
unchanged. Moreover, when it falls off as 1/a 2 , like clas- 
sical radiation or as shown to be the case for the class of 
initial state perturbations described in appendix [Dl the 
curvature perturbations due to scalar and vector pertur- 
bations fall off faster (by a factor 1/a) than those due to 
tensor perturbations, which become the dominant con- 
tribution at late times. 



VIII. GENERAL CONFORMAL FIELD 
THEORIES 

The results found for free conformal scalar fields in 
the previous sections can be straightforwardly general- 
ized to any CFT for the matter sector (even strongly 
coupled ones). This is because when the background 
gffi is Minkowski spacetime, the key ingredient for ob- 
taining the linearized stress tensor expectation value 



T. 



jhi [g^ + h]j , and hence the right-hand side of 

/ ren 

eq. ([7]). is entirely determined (up to a constant factor) 
by conformal as well as Poincare invariance. Moreover, 
it transforms in a relatively simple way under conformal 
transformations and can be easily extended to the case 
of metric perturbations around a FLRW background, as 
described in some more detail below. 

The stress tensor expectation value can be obtained 
by functionally differentiating the so-called CTP effective 
action S eS [g,g'\ 0-111]: 

T^[g]) = j^;S e &[g,g'} 



V~9 

The effective action S c s [g, g'] can be written as 

Se&[g,g'} = S div [g] - SdwW] + z[g,g'] > 



(59) 



g'=g 



(60) 

where the third term which includes nonlocal contribu- 
tions results from functionally integrating out the mat- 
ter fields. Performing those path integrals for the matter 
fields gives rise to UV divergences and one needs to in- 
troduce an appropriate regularization procedure; dimen- 
sional regularization is a good choice because it is com- 
patible with general covariance (and in a number of cases 
with conformal invariance as well). Such divergences in 
S[<7, g'\ can be absorbed by local counterterms in the bare 
gravitational action. For massless fields the cosmological 
constant and the Einstein-Hilbert term do not get renor- 
malized in dimensional regularization and only countert- 
erms quadratic in the curvature are necessary, which have 
been denoted by Sdiv[g] in eq. (|60|) . More specifically, the 
gravitational counterterms for a generic CFT on a curved 
background [9l[57| are given in dimensional regularization 
by 



Sdi V [g] = 



1 



d-A 



b J c^c^V^g^x 

+ b' J £ 4 ^d d x 



(61) 



where d denotes the spacetime dimension and £4 is the 
integrand of the Euler invariant defined in eq. (l30l) . Al- 
though the (regulated) bare £[<7,(/] is invariant under 
conformal transformations of the metric0 (see the ap- 



2 Theories whose classical action is conformally invariant only in 
four dimensions can give rise to a finite counterterm quadratic 
in the Ricci scalar, which simply implies a finite change of the 
arbitrary parameter f} in eq. 10. 
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pendix in ref. [Hj]), the counterterms in 5div[<?] are not, 
which is the origin of the trace anomaly. In fact, the 
constant parameters b and b' in eq. (foTj) , which take spe- 
cific values for each CFT, are also the coefficients of the 
Weyl-squared and £4 terms in the trace anomaly (the 
nonvanishing trace of the quantum stress tensor for con- 
formal fields). So far the statements in this paragraph 
are valid for an arbitrary metric g. If one is, however, 
interested in the linearized stress tensor for metric per- 
turbations around a background g^J one only needs to 
consider terms in S c ft\g,g'] quadratic in the metric per- 
turbations or, equivalently, focus on its second functional 
derivative, which is related to the two-point function of 
the stress tensor in the background geometry. 

Let us now show the universality of the linearized stress 
tensor expectation value in two steps. The first step is to 
establish it for metric perturbations around a Minkowski 
background. By requiring conservation of the linearized 
stress tensor together with Poincare and conformal in- 
variance, it has been shown [5^| that the form of the 
properly renormalized two-point function of the stress 
tensor in a Minkowski background for any conformal the- 
ory is essentially unique (up to a constant factor), and 
so is the renormalized vacuum expectation value of the 
stress tensor in the linearly perturbed metric. Note that 
for perturbations around a Minkowski background the 
second term on the right-hand side of eq. (|6Tj) is at least 
cubic in the metric perturbations and does not contribute 
to the renormalized expectation value of the linearized 
stress tensor operator: it, therefore, depends only on the 
constant b. 

The second step is to extend this result to the case 
of perturbations around a FLRW background via a con- 
formal transformation. As already mentioned above, the 
counterterms in eq. (|61[) are not invariant under a con- 
formal transformation of the metric and lead to the fol- 
lowing change of the effective action: 

S e s[g,g'] = S eS {g,g'] - (S div [g] - S div [g'}) 
+ (Sdiv[g] - Sdivlff']) . 

A conformal transformation leaves the two terms within 
the square brackets on the right-hand side of eq. (joTj) 
invariant in four dimensions and gives rise to terms of 
order (d — 4) otherwise, which amounts to a finite con- 
tribution in eq. (|61[) . The difference between the grav- 
itational counterterms in two conformally related ge- 
ometries is, therefore, finite. The difference between 
the Weyl-squared terms gives a term proportional to 
C ilvpa C liV( " T In a in the effective action, whose functional 
derivative corresponds to the first term inside the sec- 
ond square bracket on the right-hand side of eq. (|7|). On 
the other hand, since the £4 term in eq. (|61j) vanishes up 
to quadratic order in the metric perturbations for per- 
turbations around Minkowski space, only the £4 term 
for the conformally transformed metric g contributes at 
this order. Its functional derivative corresponds to the 
terms inside the first square bracket on the right-hand 



side of eq. ([7]), which are conserved because they re- 
sult from functionally differentiating a diffeomorphism- 
invariant integral. In addition to a finite contribution to 
the coefficient of the Weyl-squared term, which can be 
absorbed by redefining the renormalization scale /2, the 
coefficient of the squared Ricci scalar term can take arbi- 
trary finite values. Its functional derivative corresponds 
to the first term on the right-hand side of eq. ([7]) and it 
generates a term OR in the trace anomaly. 

In summary, the linearized semiclassical equation ([7]) 
will have the same form for any CFT and only the nu- 
merical coefficients in front of the two square brackets 
will change depending on the values of the parameters b' 
and b, respectively. 



IX. CONCLUSIONS 

Employing the method of order reduction, we have 
solved nonperturbatively the semiclassical Einstein equa- 
tion governing the dynamics of linear metric perturba- 
tions around de Sitter spacetime when the quantum back- 
reaction of conformal scalar fields on the mean geometry 
is included. Our exact solutions establish the stability of 
de Sitter with respect to general linear metric perturba- 
tions (of scalar, vector and tensor type) and extend some 
of the existing "no-hair" results for de Sitter in classical 
general relativity to the case in which the effects of the 
quantum vacuum polarization of conformal fields on the 
semiclassical geometry are included. Indeed, we confirm 
the late-time attractor character of de Sitter space (for 
geometrical properties within a region of fixed physical 
size) by showing that the perturbations of the Riemann 
tensor, which characterizes entirely the local geometry, 
fall off with an inverse power of the scale factor. 

Perturbative solutions of the semiclassical equation for 
tensor perturbations have recently been obtained [lTj and 
it was found that the correction to the classical solution 
can grow arbitrarily large for modes spending a long time 
inside the horizon. In contrast, our exact nonperturba- 
tive solutions exhibit oscillations with decaying ampli- 
tude inside the horizon and reveal a breakdown of per- 
turbation theory for long times inside the horizon due 
to secular terms that arise when expanding perturba- 
tively oscillatory factors with a perturbatively corrected 
frequency. In addition, we have considered the effects 
due to perturbations of the initial state of the matter 
fields. In fact, in order for the state of the fields to 
continue being a Hadamard state (with no unphysical 
excitations of arbitrarily short-wavelength modes) when 
metric perturbations are present at the initial time, it is 
in general necessary to correct the states that one would 
consider in the absence of perturbations and consider in- 
stead properly "dressed" states which are adiabatic on 
the perturbed geometry at sufficiently high order. 

As explained in sec. I Villi our results are applicable to 
general conformal field theories with arbitrarily strong 
self-interaction couplings. It would also be interesting to 
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extend the results for free scalar fields to nonvanishing 
masses and arbitrary curvature couplings, particularly to 
massless (or sufficiently light) minimally coupled fields, 
which typically give rise to larger IR effects in dc Sitter. 
We plan to return to this question in future investiga- 
tions. 

Besides semiclassical perturbations of the mean geom- 
etry, it is interesting to study the quantum fluctuations 
around it, which can be achieved by quantizing the met- 
ric perturbations around the mean geometry and dealing 
with them as an effective field theory. This was done in 
ref. (33j . where the two-point quantum correlation func- 
tion for tensor metric perturbations around de Sitter, in- 
cluding one-loop corrections from conformal fields, was 
calculated. A suitable method for selecting the adiabatic 
vacuum of the interacting theory was employed and two- 
point functions compatible with de Sitter invariance were 
obtained. In particular, all secular terms cancelled out 
despite being a perturbative calculation. The reason for 
that can be understood with a simpler but qualitatively 
similar example. Let us consider a system whose dynam- 
ics is governed by a time-independent Hamiltonian with 
a small interaction term, and let us focus on the following 
two-point correlation function: 



x B(t2)U(t2 } t 1 )A(t 1 )U(ti,to)\^> ) 



(63) 



where A and B are time-local operators in the 
Schrodinger picture. Evolving an arbitrary initial state 
|^o) f° r a long period (t\ — to) will require in general 
a nonperturbative calculation. However, for eigenstates 
of the full Hamiltonian (including the small interaction 
term) one simply gets a phase factor, U(ti,t )\^o) = 



-iE(ti-to) 



\l/o), which cancels exactly with the complex 



conjugate counterpart arising from the evolution of (^ol- 
Such an exact cancellation also implies a cancellation at 
every order in perturbation theory and guarantees the 
cancellation of any possible secular terms associated with 
the period (t± — to) which may arise in a perturbative 
calculation at finite order provided that an energy eigen- 
state of the full Hamiltonian is properly selected as the 
initial state. Nevertheless, additional secular terms may 
arise due to the time evolution Wfa, t\) in case of large 
time differences between the arguments of the correlation 
function. Thus, a perturbative calculation of the ground- 
state correlation function will be valid for small (£2 — ti) 
no matter how large (t± — to) is, but will break down 
for large (t% — t\). The situation is analogous for met- 
ric perturbations around de Sitter. In general a nonper- 
turbative calculation is needed to evolve arbitrary initial 
states other than the adiabatic vacuum for a sufficiently 
long time, but it is also required in order to calculate loop 
corrections to the two-point function with respect to the 
adiabatic vacuum which are valid for large invariant in- 
tervals (both for time-like or space-like separations). In 
this respect, the connection provided by stochastic grav- 



ity [60J between the solutions of the semiclassical Einstein 
equation and the two-point quantum correlation func- 
tions for metric pert urbations with (resummed) matter 
loop corrections [61j suggests that the nonperturbative 
semiclassical solutions found here, and the methods em- 
ployed to obtain them, could be exploited to compute 
matter loop corrections to the two-point function of the 
metric perturbations valid for arbitrarily large separa- 
tions. 
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Appendix A: Metric expansion 

Given the perturbed metric g M „ and expanding 
through quadratic order in the metric perturbation we 
have 



h = v ab h ab , 
f~9 = 1 + \h + h 2 - -h^ + 0(h 3 ) 

2 o 4 



(Al) 



where indices are raised and lowered with the unper- 
turbed metric rj^, i.e. we regard as a tensor field 
in flat space. For the Christoffel symbols we get 

= \s a liV - 1 -KS\ v + o{h?), 
s a MV = d IA hS + d v h«-ff*h MV . 

The calculation of the curvature tensors can be done 
straightforwardly and we obtain 
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R = (O^h^ - Dh) + (d v d M h + □ V - 29^^ V) 

- i (29 CT /i 1/CT - (2d T K T - d v h) + i (39 7 /i^ - 2^4) (cW 5 ) + 0(/i 3 ) . 



where □ = rf^dpdv and (9 M = r}^ v d v . 
The Weyl tensor is given by 



a/3 



-R5?ji 



75 (d-2)^ Us ^(d-l)(d-2) IXU ^ U ^ 

(A4) 



In four dimensions it obeys the identity 



(A5) 



which can be proved [62l . l63j by expanding the equality 



(A6) 



Alternatively, it can also be proved using the Gaufi- 
Bonnct theorem. 



Appendix B: Conformal transformation 



Under the conformal transformation 



(Bl) 



9iiv a Qin> 
the Christoffel symbols transform as 

f£, = + a" 1 {8%8l + 5^5; - g^g a °) d a a , (B2) 
and the curvature tensors become 

R a pj6 = R a i3 7 s - 2a~ 2 5'Y y gs]p{V y a){V a a) 

+ 4g aT 6^g s][T [a^VflVva - 2a~ 2 : (V m a)(V »] 
Rfiu = R^u - 2eT 1 V M V 1 ,a + 4a" 2 (V M a)(V„a) 
- gp,v [er 2 (V CT a)(V CT a) + a _1 n fl a] 



Appendix C: Special functions 



We define the entire function Ein(z) by 



1/0 fc=i 



(Cl) 



Its asymptotic expansion at infinity (for Sfte a < 0) is 
given by 



EiiK.wi - ■ !.h( iw] . '( - ) (/• -.. x ) . (C'ii 



where 7 is the Euler-Mascheroni constant. 

For completeness we list here some properties of the 
Bessel functions which are needed in this paper. Their 
limits are given by 



J n (x) -> 



1 



Y n (x) -> - 



r(n + l) V2 

T" r 

— cos X 

7TX 

r(n) /2\ n 



(a; ->• 0) , (C3) 



/ 2 r 7r 
J n (x) -> a/ — cos a; (2n + l) (x->oo), (C4) 

V 7TX L 4 



(HO), (C5) 



/ 2 r 7r 

F„(a;) — M / — sin 1 (2n + 1) (x -> 00) , (C6) 

V 7TX L 4 



and for their derivatives we have 
d 



i 2 R = R-6a,- 1 n 



3" > 



dx 
_d 

da; 



[x n J n {x)] = x n J n -i{x) , 
[x n Y n {x)]=x n Y n _ x {x). 



(C7) 



(B3) 



where is the covariant derivative associated with g^ v . 



Expanding the Bessel functions with respect to the or- 
der, we get 



—x 2 v Ja_ v {x) = — a; cos a; + sin ir + v 



2 sin a; + (a; cos a; — sin a;) (7 + In (2a; 2 )) 



Sm [Ein(2ix) (1 + ix) < 



0{y 2 ) 



(C8) 
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X 2 V Y'i_ v (x) = — cosx — XSVCLX + V 



2 cos x — (cos x + x sin x) (7 + In 2) 



!fte [Ein (2\x) (1 + ice) e lx ] — n (x cos x — sin x) 



(C9) 



^ 



Appendix D: Regular initial states 

In deriving the linearized semiclassical Einstein equa- 
tion 13 from the renormalized CTP effective action 
(42I l43j . a number of integration by parts with respect 
to the spacetime variables were performed in order to 
write the nonlocal term as it appears in eq. ([7]). After 
renormalization, the kernel H(x — x'; ft) in the effective 
action is a well-defined distribution provided that the 
metric perturbations fall-off sufficiently fast at spatial in- 
finity and at the asymptotic initial time (there is no such 
requirement for the asymptotic future due to the causal 

I 



I 

nature of the kernel). However, when giving initial con- 
ditions at some finite time 770 , the linearized stress tensor 
expectation value gets boundary contributions from that 
lower integration limit which diverge when the stress ten- 
sor is evaluated at 770. 

Working in Fourier space for the spatial coordinates, 
the nonlocal term in eq. (JTJ) can be written as 



3a 
ni a 2 (77) 



H(r] - rj', p; i^Apvir)' , p) dr/ , 

(Dl) 



where the Fourier transformed kernel is given by 



H(rj - rj',p;n) =cos[|p|(?7 



77')] d-lim 

e-i-0 



e(rj-7/-e) 



rj — rj 



+ 8(r)-r)') (ln(/ie) + 7) 



(D2) 



see ref. [33j for its computation and for further de- 
tails. The structure of the nonlocal term for each spa- 
tial Fourier mode is then analogous to that for the case 
of spatially isotropic and homogeneous metric perturba- 
tions studied in refs. [H, |23| . There it was shown in 
detail that boundary terms at 770 arise when writing the 
result in the same form as in eq. (|D1[) and that they give 

a contribution to (t$ (77)^} which diverges for rj = t/q. 

Furthermore, it was clarified that the reason for such di- 
vergences is the fact that although the Bunch-Davies vac- 
uum is a Hadamard state |9( with regular UV behavior in 
de Sitter spacetime, in general that is no longer the case 
in a perturbed geometry: with respect to well-behaved 
adiabatic vacua associated with this geometry it exhibits 
excitations of modes with arbitrarily short wavelengths. 

In refs. [22|, HH a simple method was employed for con- 
structing a family of properly "dressed" Gaussian initial 
states which are regular on the perturbed geometry. The 
states are prepared by evolving an asymptotic Bunch- 
Davies vacuum state from —00 to 770 in a given (nondy- 
namical) perturbed geometry which is asymptotically de 
Sitter and matches the dynamical geometry at 770. The 
metric perturbations during this preparation period can 
be fairly arbitrary, which allows the generation of a wide 
family of Gaussian states, but need to fulfill a few require- 
ments: they need to fall off sufficiently fast as 77 — > —00, 
so that the time integral in the nonlocal term converges, 
and they need to be small enough so that their contribu- 
tion to (f$(r))\ can be treated as a small perturbation. 
Moreover, the matching at 770 between the nondynami- 



I 

cal metric perturbations during the preparation period 
and the dynamical ones has to be smooth enough: up 
to fourth order, which is the maximum number of time 
derivatives that can appear in A^ v . In fact, as shown in 
detail in refs. [H, [HI, requiring a smooth matching up 
to this order is exactly equivalent to demanding that the 
states generated in this way are of fourth adiabatic order, 
the standard requirement for regular states with a finite 
renormalized stress tensor expectation value Q. 

When the preparation method described in the pre- 
vious paragraph is employed, the nonlocal contribution 
in eq. (|D1[) can be naturally separated into two contri- 
butions which result from the following splitting of the 
time integral: 



d7;'= / d?/+ / drj'. 

-OO J 770 J —OO 



(D3) 



The first term on the right-hand side will contain the dy- 
namical metric perturbations in the integrand and it will 
vanish when using order reduction since A^ vanishes in 
that CcLSG, cLS described in sec. IIV CI On the other hand, 
the second term will give a finite contribution in eq. (ID1[) , 
even for 77 = 770, provided that A^ is regular and vanishes 
at 7; = 770 , which follows from the condition of sufficiently 
smooth matching required above. This term gives a con- 
tribution to the right-hand side of eq. ([7]) which can be 
interpreted as a perturbation of the stress tensor expec- 
tation value associated with the modified initial state. In 
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fact, the equation can then be rewritten as eq. (|47|) with 

ST^ri, p) = ^-y J Ap V {ri , p) H(r) - 77 , p; /x) cfy . 

(D4) 

The integral is finite for 77 > 770 (remember that the met- 
ric perturbations are required to fall off sufficiently fast 
as j] — \ — 00 so that the lower integration limit is conver- 
gent) and so is its limit 77 — > 0. Therefore, at late times 
ST^v decays like 1/a 2 , clearly fulfilling the requirement 
of sees. [yn^Hyn Efl so that the inhomogeneous solutions 
of the semiclassical equation do not alter the conclusions 
about the semiclassical stability of de Sitter spacetime in 
this context. 



and (|21[) . have not appeared explicitly in the literature 
before. On the other hand, eq. (|2"21 for the tensor per- 
turbations can be directly compared with Starobinsky's 
eq. (7) in ref. [35]. To do so, we need to calculate first 
the curvature terms for the unperturbed physical metric 
9^v = a 2 Vni> which were employed by Starobinsky, and 
which are given by 

Bl = 3H 2 , R=12H 2 , R'=0. (El) 

If we insert the following result for the kernel in Fourier 
space obtained in ref. [42j: 



Appendix E: Comparison with Starobinsky's 
equation 

To our knowledge, the semiclassical equations of mo- 
tion for the scalar and vector perturbations, eqs. ([20)1 



H(x;n) = -5 



In 



4 -i7re(V)sgnp le i ^ 7 ^, 

(E2) 

employ the mode decomposition introduced in eq. (|33j) 
and divide by a 2 , eq. (|2"2")l becomes 



[l - an 2 H 2 - 2(3k 2 H 2 ] (c£ + k 2 g±) + 2 [l - an 2 H 2 - 2(3k 2 H 2 ] Hag' ± = 
2Ha(g^ + k 2 g' ± )+H 2 a 2 (g'i-k 2 g ± ) 



3an 2 a 2 



(E3) 



-ik u r, 



g ± (v')e ik "^ dry' 



2 ,2 



(k 2 ) 



hi 



fi 2 a 2 



+ me(~k 2 )sgnk° 



dk° 
~2^ 



This coincides exactly with Starobinsky's equation if 
we take into account that his constants M 2 , H 2 and 
G£/(60n) correspond to 2/(/3k 2 ), -2/(aK 2 ) and 3aK 2 , 
respectively. Note that Starobinsky's H 2 is positive be- 
cause he considers photon fields, while we have conformal 
scalars, for which the corresponding constant changes 
sign |64| . 



In Starobinsky's inflationary model the expansion is 
driven by the trace anomaly of the quantized matter 
fields, so that the actual value of Starobinsky's Hubble 
parameter H is different from our H , where the expansion 
is driven by the cosmological constant A. Nevertheless, 
this has no effect on the form of eq. (|E3[) . 
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